Search for: All records

We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo by applying automatic differentiation (AD). Biases and statistical efficiency of the resulting estimators are discussed. Our approach demonstrates that AD enables the calculation of reduced density matrices with the same computational cost scaling per sample as energy calculations, accompanied by a cost prefactor of less than four in our numerical calculations. We investigate the role of self-consistency and trial orbital choice in property calculations. We find that orbitals obtained using density functional theory perform well for the dipole moments of selected molecules compared to those optimized self-consistently.

 

We present efficient algorithms for using selected configuration interaction (sCI) trial wave functions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC). These advances, geared toward optimizing computational performance for longer configuration interaction expansions, allow us to use up to a million configurations in the trial state for ph-AFQMC. In one example, we found the cost of ph-AFQMC per sample to increase only by a factor of about 3 for a calculation with 10 4 configurations compared to that with a single one, demonstrating the tiny computational overhead due to a longer expansion. This favorable scaling allows us to study the systematic convergence of the phaseless bias in auxiliary field quantum Monte Carlo calculations with an increasing number of configurations and provides a means to gauge the accuracy of ph-AFQMC with other trial states. We also show how the scalability issues of sCI trial states for large system sizes could be mitigated by restricting them to a moderately sized orbital active space and leveraging the near-cancellation of out of active space phaseless errors.